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FOREWORD 


A  major  problem  in  all  engineering  design  is  that  of  reliability  and 
quality  control.  Standard  procedures  exist  for  assuring  the  desired  reli¬ 
ability  under  ordinary  circumstances.  This  report  describes  an  unusual 
situation,  in  which  the  tolerances  on  the  components  of  an  assembly  are 
not  small  enough  to  ensure  that  the  assembly  wild,  work  properly,  and  provides 
estimates  of  the  probability  of  a  defective  assembly  in  this  case. 
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ABSTRACT 


This  report  describes  an  investigation  of  how  errors  in  components  of  an 
assembly  can  affect  its  performance.  In  particular  the  report  de:  with  the 

situation,  uncommon  in  engineering  practice,  where  the  output  tolerance  of 
the  assembly  may  be  violated  even  tliougn  the  tolerances  on  the  components  are 
all  met.  This  situation  is  analyzed  to  estimate  the  probability  that  the 
output  tolerance  will  be  satisfied  given  that  the  component  tolerances  are 
met.  Three  methods  are  described  for  estimating  this  probability,  their 
results  eire  compared  in  a  number  of  cases,  and  a  best  method  is  chosen. 
Several  simple  rules,  suitable  for  preliminary  estimates,  are  also  given. 

An  example  is  worked  out  showing  a  simple  application  of  the  method. 


1.  INTRODUCTION 


This  report  deals  with  certain  aspects  of  the  general  problem  of  errors 
and  tolerances  in  the  design  and  testing  of  equipment.  It  is  presumed  that 
the  piece  of  equipment  is  required  to  operate  at  a  certain  level  of  output. 
Ordinarily  the  designer  assigns  a  certain  error -tolerance  to  this  output, 
chosen  so  that  the  equipment  will  function  properly  if  the  output  error 
satisfies  its  tolerance.  The  output  error  usually  arises  from  errors  in  the 
various  components  that  have  been  assembled  to  make  the  piece  of  equipment. 

The  designer  will  customarily  know  the  relation  between  the  output  error  and 
the  component  errors.  Common  practice  (see  Bowker  and  Liebermarf ^ )  is 
that  the  designer  will  combine  this  relation  with  the  output  tolerance  to 
find  tolerances  on  each  component  such  that  satisfaction  of  these  component 
tolerances  will  ensure  that  the  output  tolerance  is  met. 

We  are  concerned  here  with  the  uncommon  situation  where  satisfaction 
of  the  component  tolerances  does  not  ensure  satisfaction  of  the  output 
tolerance.  This  state  of  affairs  can  arise  when  an  error  has  been  made  in 
choosing  the  component  tolerances,  or  when  it  is  impractical  (or  too  expensive) 
to  make  the  component  tolerances  small  enough.  In  either  case  we  must  face 
the  possibility  that  all  the  components  will  meet  their  tolerances  but  some 
of  the  assembled  pieces  of  equipment  will  not  work  properly.  The  practical 
information  that  we  want  is  the  probability  that  the  output  tolerance  will 


be  satisfied.  With  this  information  we  can  estimate  how  many  extra  pieces 
of  equipment  must  be  manufactured  on  the  average  in  order  to  obtain  a  given 
number  of  workable  assemblies. 

In  the  following  section  we  shall  describe  the  general  procedure  for 
estimating  the  probability  that  the  output'  satisfies  its  tolerance,  supposing 
that  each  component  error  is  normally  distributed  with  zero  mean,  known 
variance  and  known  tolerance,  Three  mathematical  methods  are  given  for 
carrying  out  the  calculations.  One  is  of  Monte-Carlo  type  and  is  described 
in  Section  3*  The  other  two  methods  use  the  Characteristic  Function  in 
different  ways.  Section  4  gives  formulas  for  the  Characteristic  Functions  of 
the  various  distributions,  and  Sections  5  sad  6  use  these  formulas  in 
estimating  the  desired  probability.  Various  simple  approximations  and  limit¬ 
ing  cases  are  examined  in  Section  7*  Section  8  describes  the  results,  which 
are  then  discussed  in  Section  9.  The  report  closes  with  a  simple  example  of 
how  these  estimates  might  be  used  in  practice.  Section  10. 
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We  let  Yj  be  the  error  in  the  j-th  component,  js  1, , 
and  XQ  is  the  error  in  the  output.  The  relation  between  the  output  error 
and  the  component  errors  is  taken  as  linear, 

N 

Xo= 

J=i  J 

where  the  Cj  are  assumed  to  be  known  constants. 

It  is  assumed  initially  that  Yj  is  normally  distributed  with  zero  mean 
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and  variance  Oj  .  Then  we  may  define 


Xj  =  CjYj 
Sj-lCjlcj 


j=  I,  z.,...  w 


and  the  relation  (l)  can  be  written 

M 


XQ=  Z  Xj 

J=l  J 


where  Xj  is  normally  distributed  with  zero  mean  and  variance  Sj  • 

We  let  Dj  >0  be  the  tolerance  on  the  error,  Yj,  in  the  j-th  component, 
and  B0^  0  be  the  tolerance  on  the  output  error,  XQ.  Thus,  when  Yj  satisfies 


its  tolerance, we  have 


IVjNDj 


(1) 


(2) 


(3) 


3 


and,  if  X0  satisfied  its  tolerance,  then 


IX.NR.  - 

We  define  also 

as  the  tolerance  on  Xj,  so  that,  if  Xj  satisfied  its  tolerance,  then 

We  notice  also  that 

Bj/Sj  =  Dj/Oj 

We  now  define  certain  probabilities.  Pj  is  the  probability  that  the 
j-th  error  satisfies  its  tolerance,  i.e., 

Pj  -  prob^N  Dj]=  Prob[lX^Bjl 

Pc  is  the  probability  that  all  ccmponent  errors  satisfy  their  tolerances.  We 
assume  that  the  component  errors  are  independent  of  each  other,  and  therefore 

prok^lXjl^Bj  for  all  j] 

N 


Further,  ve  define 


PQ  =  prob  £~- the  output  error  satisfies  its  tolerance  and  all  component 
errors  satisfy  their  tolerance sj 
P0  =  prob  £“  |X0|6  60  and  jY^l  ^  0>j  for  all  j  J 
The  theorem  on  compound  probability  asserts  that 

P0  =  prob  c  w  6  &„ ,  given  that  IZ  I  ^  Bj  for  all  j  J 

x  prob  Z”  Bj  for  all  j  J 

The  probability  that  is  of  greatest  practical  interest  is 

P*  =  prob  f  \X0\^  B0  ,  given  that  >*iU  Bj  for  all  j J 

Then  we  can  write  (9)  with  the  aid  of  (7)  and  (10)  as 

p*=  P./Pe  =  p„/|r't>- 

Finally,  it  is  useful  to  define 

A  =  I  B,  =  llCjlDj 
j=i  J  J=i 

If  all  the  components  satisfy  their  tolerances,  then,  using  a  familial* 
property  of  inequalities,  we  find 

J 


(10)  j 


and  so,  because  of  (12),  X0  must  satisfy  the  inequality 


tXUA 
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If  A  4  6,  then  (13)  implies 


IX0UA«  Bo 


In  this  case  we  see  that,  if  each  component  satisfies  its  tolerance,  the 
output  error  X0  must  always  satisfy  its  tolerance,  and  from  (10)  and  (ll) 
we  conclude  that 


This  case  is  the  common  one  in  design  practice,  i.e.,  the  tolerances  are  set 
so  that,  if  each  component  meets  its  tolerance,  the  output  will  necessarily 
satisfy  its  tolerance.  However,  in  this  paper  we  are  interested  in  the 
opposite  case,  where 

A>  60 

and 

o-<  1 


Our  main  objective  is  to  estimate  P*.  We  define  F*  (XD)  as  the  density 
function  of  the  output  when  the  separate  component  errors  all  satisfy  their 
tolerances.  Since  the  component  errors  are  normally  distributed,  their  density 
functions,  when  they  satisfy  their  tolerances,  are  symmetrically-truncated 
normal  distributions,  and  F*  (Xq)  is  the  density  function  of  a  finite  sum  of 
such  distributions.  P*  is  the  integral  between  -B0  and  B0  of  F*  (X0). 
Unfortunately  F*  (Xq)  is  not  easily  expressible  in  terms  of  the  parameters 
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of  the  component  density  functions.  However,  we  can  make  a  number  of  simple 
comments  about  tne  behavior  of  F*  (Xo). 

(i)  If  the  component  error  tolerances  are  all  very  large,  i.e., 
e>j»  Sj  ,  each  component  error  is  approximately  normally  distributed, 
hence  F*  (X0)  is  approximately  a  normal  function. 

(ii)  If  N,  the  number  of  error  components,  is  large,  the  Central 

Limit  Theorem  leads  us  to  expect  that  F*  will  be  approximately  a  normal  function. 

(iii)  Contrariwise,  F*  will  depart  furthest  from  normality  when  some 
Bi/Sj  are  small  and  when  N  is  small.  This  will  be  particularly  so  when  one 
component  dominates  all  the  rest,  so  that  effectively  N  =  1. 
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3.  MONTS -CARLO  METHOD 


This  method  is  based  on  counting  the  numbers  of  successes  in  sampling 
from  distributions  that  are  random,  independent  and  normally  distributed 
with  zero  mean.  A  computer  program  was  written  to  carry  out  this  procedure. 


The  program  has  as  input  the  quantities  Cj,  Dj  (j  =  1,  2,  ...  N)  and 
B0,  together  with  the  list  of  random  numbers.  L  sets  of  N  random  numbers  are 
read  in  successively.  For  each  set,  the  N  random  numbers  are  taken  as  the 
values  of  Yj  (j  =  1,  2,  ...  N),  and  the  value  of  XQ  is  calculated  from  (l). 
Counts  sure  made  of  the  following  quantities: 


nJ  ^ 

nc  is 

nQ  is 


=  1,  ...  N)  is  the  number  of  cases  for  which 

the  number  of  cases  for  which 

the  number  of  cases  for  which  IXil*  fc 

and  !X0|£$0 


'XjUDj 


for  all  j. 

for  all  j 


Then  we  obtain  the  estimates 
Pj  =  nj/L 


Pc  ~  nc A* 


Pq  -  no /L 


P*  -  Pq/Pc  -  nQ/nc. 


j  =  1,  2,  ... 


N 


In.  using  this  method  we  had  to  choose  L  large  enough  so  that  reasonably 
stable  estimates  of  P*  were  obtained.  The  choice  L  =  200  was  used,  and  the 
entire  procedure  repeated  four  times  with  different  sets  of  random  numbers. 

The  final  estimates  of  the  probabilities  are  given  as  the  means  of  the  results 
for  the  four  repetitions. 
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4.  CHARACTERISTIC  FUNCTIONS 

The  remaining  two  methods  of  estimating  P*  employ  the  Characteristic 
Function  (or  Fourier  Transform)  as  the  main  tool  in  the  analysis.  In  this 
section  we  present  the  general,  formulas  that  form  the  basis  of  t^ise  methods. 

If  F  (x)  is  a  density  function,  then 

Xa+«? 

erAty  foo  <JX 


<^>(0  = 


Xs -QE? 


(14) 


is  its  Characteristic  Function  or  Fourier  Transform.  It  is  unnecessary  to 
dwell  on  the  properties  of  the  Characteristic  Function  which  are  well-known. 

We  record  only  one  formula,  which  is  easily  derived  from  the  Complex  Inversion 


Relation, 


CO 


j  Fto)(U  =  Tr'1  j  sinftiQAX 


-t=-oo 


This  formula  expresses  the  area  under  the  density  curve  between  -X  and  X 
(or  the  probability  that  the  variable  lies  between  -X  and  X)  directly  as  an 
integral  of 

We  define  Fj  (Xj)  as  the  density  function  for  Xj  and  as  the 

corresponding  Characteristic  Function.  Similarly  is  the  Character¬ 

istic  Function  of  the  output  distribution,  F*  (X0).  Since  the  Xj  are  assumed 
to  be  independent,  it  is  well-known  that 
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(15) 


(16) 


Both  methods  of  estimating  P*  are  based  on  finding 
by  means  of  (l6). 


tw  from  the 


To 


find  j  (t;)  we  first  write  down  the  density  function  of  Xj, 


which  is  that  for  a  normal  distribution  with  variance  Sj  ,  truncated  at 


±Bj> 


j [V/(^  del 


=  0 


This  is  then  substituted  into  (l4),  and  we  find,  after  some  manipulation,*- 

w  =  ( ri + xpS)} 

4 


where 


-- 


6jA5j-\ff)  >o 

Pi  =  tSj/VT 

A  series  representation  of  V+)  in  real  terms  may  be  derived  by 
expanding  the  Error  Function  in  a  Taylor  Series  about  .  By  means  of 
(A.  5)  we  find 

1  %  4  jJ  br01  J 


(17) 

(13) 

(19) 


(20) 


*For  completeness  a  list  of  the  basic  forma  as  relating  to  the  Error 
Function  is  given  in  Appendix  A. 


11 


■where  Hk  is  the  Hermite  Polynomial  of  degree  k.  This  series  converges 
absolutely  for  any  finite  values  of  /»J  and  V*  .  Further,  we  may  expand 
6.  about  =  0  and  obtain  explicitly  the  leading  terms  (up  to  t^) 

in  the  expansion  of  fjW  about  t  =  0, 

ftw- 1  -  V(1  -  *C*f%v*0 

+  s;{3~  1tt-V2 


Although  the  series  in  (20)  converges  for  any  finite  values  of  and 
Yj  ,  the  convergence  is  slow  when  is  large,  and  an  alternative  method  of 
computation  is  needed.  For  this  purpose  it  is  convenient  to  use  the  real 
and  imaginary  parts,  Wr  and  W±,  of  the  complex  function,  W,  defined  in  (A.9). 
From  (17),  (A.  10)  and  (A.ll)  we  obtain  the  exact  formula 

e7?  £  { ^  " e" ^  '-V' 

A  rational  af  roximation  for  W  is  given  in  (A. 12),  and  from  it  we  may  derive 
the  following  approximations  for  Wr  and  W^: 

■vUfi  I  rK  {(-  y.  +/>. 

K=\ 


(24) 


Here 


» 


and  X~K  and  Tj^  are  constants  of  the  approximation,  listed  in  Appendix  A. 
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5.  APPROXIMATION  USING  MOMENTS 


This  method  consists  of  assuming  that  the  output  density  function  is 
approximately  of  the  form 

g  (X.)  -  S' WT.V'Wr  +  Q^XyS?]  e  X°Va5l)  (26) 

2 

vhere  G0,  Gg  and  G4  are  constants  to  he  determined,  and  S  is  the  exact  second 
moment  of  the  output.  The  constants  G0>  G2  and  G]^  are  chosen  by  matching 
moments,  i.e.,  by  using 
00 

^  X0ZK^Uo)clX0  =  M1K  ,  0,lfz  (27) 

-00 

where  %  are  the  exact,  even-ordered  moments  of  the  output  distribution. 

M0,  Mg  and  M^  are  determined  by  using  the  well-known  relation 

= ',-ir  4^  c°>  <*> 

Because  the  If)  are  characteristic  functions  and  are  all  even  in 

t,  we  have 


4^(0)  =  1 ,  =  4>j'"(o)  3=0 


Differentiating  (l6)  and  using  (29)  we  find 


=  4>„(o)  -  1  •  (30) 

(31) 

~  b"(°)  =  Z  [-  VloSj  (32) 

^3  '  -*•  'f',,  (o)  =  0  (33) 

H*=  3[t,/,fo)]  +  Z[  t)B(o)  “  *M]  (34) 

From  (2l)  we  obtain 

-  *>>-  ^{i-L^A'(zyA(^)]}  (35) 

^°)=  ${'  tij1)}  (36) 

where 

Aj  =  Bj/Sj  =  2!  it]  (37) 


and  the  functions  A  and  A'  are  defined  in  (A. 2)  and  (A. 6).  We  calculate 
M2  and  Mi|  by  inserting  these  values  in  (32)  and  (34). 


If  we  combine  (26)  with  (27)  and  use  the  general  formula 

00 


C2-rr)V1jvlneV/<1S^v=  S^tVcik-I), 


—CD 


(38) 


we  obtain  the  following  linear,  algebraic  equations  for 
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G0,  G2  and  G^: 


+  =  1 

Q„  +  3^+ls-&H  = 

3&„  +  \os-^=  M4/Sm=  Mh/*C 


This  system  is  solved  for  G0,  G2  and  G^,  the  results  are  inserted  into 
(26),  and  we  can  then  calculate  the  approximate  P*  by  means  of 

P'=  S  »  C0iXo=  l\  j  (X„)dX0 


Using  the  general  relation  (38)  again  we  find  the  following  formula  for  P* 


where 


p*(2)  =  A©  -  HI  (3  - 11)  A' (2) 

(40) 

H  -  {3-(MH/Mx)}Alt 

(41) 

2  =  6„/S 

(42) 

We  shall  call  the  estimate  of  F*  given  by  (40)  the  moment  approximation. 
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6.  NUMERICAL  INTEGRATION  OF  THE  CHARACTERISTIC  FUNCTION 


i.e. 


This  procedure  consists  merely  of  carrying  out  the  integration  of  (15), 
,  evaluating 


co 

P*  =  iF\x.A=(2/lr^t(t)-r1SvnC6t>t)  it 


-6. 


where  is  calculated  from  4,(0  by  (l6).  The  4j(« 
are  evaluated  by  use  of  (20)  vhen is  of  moderate  size  and  (22)  -  (25) 
when  p j  is  large. 


(43) 


In  general  it  is  necessary  to  evaluate  the  integral  by  numerical  means. 
To  carry  this  out  with  sufficient  accuracy  is  sometimes  difficult  because 
tw  dies  away  in  oscillatory  fashion  as  "t— *<JD  .  Usually,  most  of 
the  contribution  to  the  value  of  the  integral  comes  from  near  t  =  0,  but 
significant  contributions  can  also  come  from  further  out,  where  <Ut) 
may  oscillate  rapidly.  In  cases  where  a  significant  contribution  comes  from 
the  region  of  rapid  oscillation,  the  integration  must  be  carried  out  with 
great  care. 


The  following  procedure  was  adopted  for  carrying  out  these  integrals. 
A  fundamental  range  of  t  say  T,  was  chosen,  roughly  small  enough  so  that 


has  no  more  than  5  zeroes  in  the  range  (n-l)T  itinT  for  n  =  1,  2,  3,  ..  M. 
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M  is  taken  so  large  that  the  ranges  beyond  n  =  M  contribute  negligibly  to 
the  integral.  For  each  range  of  t  the  integral  was  evaluated  by  Gaussian 
Integration,  and  the  total,  integral  obtained  by  adding  the  results  for  ail 
the  ranges.  Seme  experimentation  was  needed  to  find  a  suitable  number  of 
points  to  use  in  the  Gaussian  Integration  and  to  determine  how  large  M  should 
be. 


The  computer  program  that  carried  out  this  evaluation  of  P*  was 
occasionally  slow-running  and  was  used  primarily  for  spot-checking  the  results 


of  the  other  methods 


7.  SIMPLE  APPROXIMATIONS  AND  LIMITING  CASES 


A  number  of  obvious,  simple  approximations  for  P*  can  be  derived,  and 
we  describe  three  of  them  briefly  here. 

(i)  If  all  the  component  tolerances  are  very  large,  i.e.,  if 


VL  *  Vs,  »  1  + 


or  <»•"  j 


the  distribution  associated  with  each  Xj  is  approximately  normal  with  variance 
.  Then  F*  (XQ)  is  approximately  a  normal  distribution  with  variance 


s  x 

,  where 


|si 


S1  =  7  Sj 


(bk) 


Hence  we  may  write  this  estimate  of  P*  as  a  function  of  B0  in  the  form 


ft  (VSt)=  A  (BA) 

(ii)  A  different  approximation  may  be  obtained  if  we  assume  that 
F*  (XQ)  is  a  normal  distribution  with  variance  ,  truncated  at  points 


i.e. 


F*(x.)  =  o 


\K\  5  T, 


(45) 


In  this  case  we  find  the  following  estimate 


P*  (6./^  =  A  ( tp/Sjj)/ A  (r.  /  5,) ,  B.<  tc 


=  1 
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1  i 


The  properties  of  the  error  function  imply  that 


The  accuracy  of  this  approximation  depends  on  how  well  F*  (XQ)  is 

approximated  by  a  truncated  normal  distribution  and  how  precisely  we  can 

estimate  There  are  two  cases  in  which  this  approximation  may  be 

tolerably  accurate.  First,  if  one  component,  say  the  L-th,  dominates  the 
£  2L.  X  *  * 

others,  i.e.,  »  the  output  will  be  approximately 

that  of  the  dominant  component,  X^,  which  is  a  normal  distribution  truncated 
at  +1}^.  In  this  case  we  expect  that 


with  TQ  =  B^,  will  be  a 


fair  approximation  to  P*.  Second,  if  all  the  C.  are  roughly  equal,  and  all 

J 

the  Zj  are  roughly  equal  to  Z^,  say,  then  we  expect  that  F*  (XQ)  will  be 
approximately  a  normal  distribution  with  variance  ^  truncated  at 

= 

Hence  in  this  case  also  ^<6A)  should  be  a  decent  approximation 

to  P*. 


(iii)  A  third  simple  approximation  is  obtained  by  setting  G^  =  0  in 

(26)  and  choosing  GQ  and  Gg  such  that 

OO 

j  =  M,k 

Then  we  get  as  an  approximation  for  P*  merely  the  first  term  of  (42), 

?*&)=  A®. 


4t>r-  K  =  Oj  1 


(46) 


(47) 
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8.  RESULTS 


All  the  cases  discussed  here  have  four  components,  i.e.,  N  =  4,  and  all 
have  (IT  *  1,  j  =  1,  2,  3>  4.  Eleven  different  cases  were  studied  with  various 
values  for  the  Di  and  C  .•  as  shown  in  Table  1.  For  each  case  results  were 

J  J 

obtained  in  the  form  of  graphs  of  P*  as  a  function  of  B0/Sj  and  are  displayed 
in  Figures  1-12..  Each  graph  shows  the  mean  and  standard  deviation  of  the 
four  repetitions  of  the  Monte -Carlo  Method  as  well  as  the  moment  approxima- 
tion  for  that  case.  The  values  of  Sj,  Mg  =  S  and  M^  are  listed  in  Table  2. 

The  method  of  integrating  the  Characteristic  Function  was  used  only  at 
points  where  sizeable  discrepancies  were  found  between  the  Monte-Carlo  x\'sults 
and  the  moment  approximation.  These  points'  are  shown  on  the  appropriate 
graphs  and  compared  with  the  other  methods  in  Table  3» 

Several  additional  graphs  shew  comparisons  among  the  moment  method 
predictions  for  different  cases.  The  comparison  among  Cases  (I),  (II)  and 
(ill)  is  shewn  in  Figure  12,  Cases  (Vi),  (VIl)  and  (VIII)  in  Figure  13  and 
Cases  (IX),  (X)  and  (Xl)  in  Figure  14.  Also,  the  simple  approximation 
P-j-*  is  shown  in  Fig-are  12,  and  PII*  in  Figures  13  and  l4  for  relevant  values 
of  T  . 


9.  DISCUSSION 


We  shall  first  compare  and  comment  on  the  results  obtained  by  the  various 
methods,  then  suggest  procedures  for  estimating  P*  under  various  circumstances 

Figures  1  to  11  show  that  the  agreement  between  the  Monte -Carlo  method 
and  the  moment -approximation  is  reasonably'  good  in  a  general  sense.  We  see 
from  Table  3  that,  when  the  results  do  not  agree  well,  the  integration  of 
the  characteristic  function  almost  always  agrees  with  the  moment  approximation 

Of  the  three  methods  one  expects  that  integrating  the  characteristic 
function  should  be  the  most  accurate.  The  Monte -Carlo  method  is  usually 
thought  to  be  somewhat  inaccurate  unless  a  very  large  number  of  samples  is 
used,  and  the  above  results  suggest  that  this  is  the  case  here.  Procedures 
like  the  moment -approximation  are  fairly  common  in  statistics  and  often 
give  satisfactory  accuracy.  However,  there  are  two  theoretical  defects  of 
the  moment-approximation  here  that  are  worth  mentioning.  First,  the  approxi¬ 
mate  density  function  g  (XQ)  is  continuous  (see  Equation  (26))  but  the  true 
density  function,  F*  (X0),  is  discontinuous.  In  fact 

F*(X0)ho 

when  X„>&,  .  Second,  g  (X0)  is  slightly  negative  for  X0  sufficiently 

large  in  many  cases.  Neither  of  these  defects  seems  to  cause  serious  errors 
in  the  estimate  of  P*  for  the  casea  studied  here  since  the  results  agree 
well  with  the  integration  of  the  characteristic  function.  If  serious  errors 
sire  to  arise,  one  would  expect,  to  find  them  when  there  is  a  single,  dominant 
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component  with  a  low  tolerance  on  it,  as  in  Case  (VIIl).  Table  3  confirms 
this  expectation,  for  we  see  that,  when  B0/Sj  =  .9^5*  the  characteristic 
function  and  moment-approximation  differ  by  .011,  whereas  the  worst  error 
observed  in  the  other  cases  of  Table  3  is  only  .001.  However,  even  in  this 
most  unfavorable  case  the  error  in  the  moment-approximation  is  small  enough 
to  be  unimportant. 

Figure  12  shows  how  curves  for  P*  change  as  the  common  tolerance  value 
for  the  four  components  increases  from  Bj/Sj  =  1  through  1.5  to  2.  As  we 
expect,  the  curves  become  lower  and  tend  toward  the  normal  curve,  given  by 
Pj*,  with  increasing  component  tolerances. 

The  effect  on  P*  of  an  increasingly  dominant  component  is  displayed  in 
Figures  13  and  l4.  Figure  13  shows  the  case  where  the  increasingly  dominant 
component  has  a  smaller  tolerance  than  the  other  components.  As  increases 
from  1  through  2  to  5j  the  curve  of  P*  is  raised  toward  the  curve  for  Jjj* 
truncated  at  B0/Sj  =  1,  to  which  it  must  ultimately  tend.  In  contrast 
Figure  l4  shows  what  happens  when  the  increasingly  dominant  component  has  a 
higher  tolerance  than  the  others.  As  increases  from  1  through  2  to  5, 
the  curve  for  P*  is  lowered  toward  the  curve  for  Fjj*  truncated  at 
B0/Sj  =  2,  to  which  it  ultimately  tends. 

When  all  the  C  j  are  roughly  the  same,  we  may  also  inquire  about  the 
effect  of  changing  the  component  tolerances  but  keeping  the  average  component 
tolerance  constant.  Comparing  Cases  (ll) ,  IV)  and  (v)  we  c^e  that  this  has 
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scarcely  any  effect  on  P*.  In  other  words,  when  the  C  j  are  roughly  equal, 
the  mean  component  tolerance  has  a  considerable  influence  on  P*  (see  Figure 
12)  but  the  variance  in  the  component  tolerances  has  negligible  effect. 

A  reasonably  extensive  comparison  of  the  simple  approximations,  (47), 
(48)  and  (49)  with  the  more  accurate  calculations  suggests  the  following  as 
a  rule: 

(i)  Use  Pj-r*  if  one  component  dominates  greatly. 

(ii)  Use  Pj*  if  Bj/Sj  >  2  for  all  j. 

(iii)  Use  Pjjj*  if  no  one  component  dominates. 

Use  of  this  rule  will  give  fair  results,  perhaps  suitable  for  an  initial 
estimate.  The  simplest  accurate  procedure  is  the  moment-approximation, 
given  by  (42)  -  (44). 


10.  WORKED  EXAMPLE 


* 


A  certain  piece  of  mechanical  equipment  is  supposed  to  operate  at  a 
load  of  900  lbs.  Three  components,  standard  items,  are  assembled  to  form 
this  mechanism.  Components  1  and  3  are  springs  and  component  2  is  an 
electrical  switch.  Their  component  errors,  Y]_,  Y2  and.  Y^,  are  related  to 
the  error  in  the  output  load  by  (l)  where 

Cj_  =  ?.40  lb/inch 
C2  =  15  lbs/volt 
n  210  lb /inch 

From  information  about  the  manufacture  of  the  components  we  know  that  the 
distribution  of  their  errors  is  roughly  normal  with  zero  mean  and 

#1:  60$  are  acceptable  at  a  tolerance  of  . 1-inch 

#2:  70$  are  acceptable  at  a  tolerance  of  1  volt 

#3=  92$  are  acceptable  at  a  tolerance  of  ,1-inch 

From  Tables  of  the  normal  distribution  we  find  the  standard  deviation  CTJ 
as  follows: 

#1:  .8  acceptable  corresponds  to  .1  =  1.28  (7^ 

(Tp  =  .O78I"  inches 

#2:  .7  acceptable  corresponds  to  1  =  1.04  0^ 

0^  =  .962  volts 
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#3:  .92  acceptable  corresponds  to  .1  =  1.75 

=  .0571  Inches 

Then 

Sx  =  240  x  .0781  =  18.7  lbs. 

52  =  15  x  *9^2  =  i4.4  lbs. 

53  =  210  x  .0571  =  12.0  lbs. 

The  tolerances  established  on  the  components  are 

=  .15  inches,  Dg  =  2  volts,  =  .07” inches 

and  therefore  from  (4)  and  (12) 

Bi  =  36  lbs,  B2  =  30  lbs,  B3  =  14.7  lbs. 

A  =  80.7  lbs. 

The  tolerance  on  the  load  needed  to  activate  the  mechanism  is 
B0  =  30  lbs. 

i.e.,  we  know  it  is  possible  that 
but  the  tolerance  on  the  load  will 

be  violated. 

First  we  find  a  quick,  rough  estimate  of  P*.  In  order  to  determine 
which  one  of  the  simple  estimates,  (45)  -  (47),  is  best,  we  find  from  (37) 

6,/S,  =\.11  ,  2.01,  \.xx 


Since  B0  <A  we  know  that 

each  component  will  satisfy  its  tolerance 
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These  are  not  all  2.  Also  none  of  t.._  Cj  dominates  all  the  others.  The 
rale  stated  at  the  end  of  section  9  suggests,  therefore,  that  we  use  Pjjj* 

To  find  Pm  we  need  to  find  M2  by  means  of  (32)  and  (35)*  Prom  the  tables 
of  the  normal  function  we  find 

A'(i,)=.uo,  A'(^=.on,  A'(£3)=  ,3T} 

Putting  these  values  into  (35);  then  combining  the  results  with  (32)  we  get 

M,=  4$u  '  s  -  U-UI 

From  (42)  Z  =  30/2I.9  =  1.37*  Using  the  estimate  (47),  the  normal  function 
table  gives 

PjH*  =  A  (1.37)  =  .829 

The  more  accurate  approximation  (4o)  involves  finding  in  addition  to 
Mg.  The  calculation  of  by  (34)  and  (36)  leads  to  =  625,000.  Putting 
this  into  (4l)  and  combining  with  (40)  we  find 

P*  (z)  =  .829  -  .006  =  .823 

Thus,  if  all  components  satisfy  their  tolerances,  the  probability  that 
the  mechanism  will  operate  at  a  load  between  870  and  930  lbs.  is  about  .82. 
On  the  average,  therefore,  if  we  need  1,000  workable  mechanisms,  we  should 
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APPENDIX  A  -  FORMUIAS  RELATING  TO  THE  ERROR  FUNCTION 


The  following  are  fundamental  formulas  relating  to  the  Error  Function 

2 

and  are  taken  from  Reference  . 

2  „ 

erf  (Z)  =  J  e"  u  Au  <a.i) 

O 

In  (A.l)  the  integration  may  be  carried  out  along  any  path  in  the  complex 
U -plane  connecting  U**0  and  U=Z.  An  alternative  definition  is 

A(Z)=  erf  (2* 'A2L)  «  (2-Tt)“'A  \  e~Vl  (*.2> 

X=-i 

We  have  also 

erf(-t)  =  -  erf  It)  (a.3) 

erf  (5)  -  erf  (¥)  (a.4) 

where  the  bar  denotes  uhe  complex  conjugate. 


Various  derivatives  of  these  quantities  can  be  found  from 

.  K+\  U  .  _ ~7^- 


(ev4  £)=  (rl)K  2rn-“V/ze  * 


K 


A'(*)=  M/il  =  (lAr)’/ze''iyi 


where  (Z)  is  the  Hermite  Polynomial  of  k-th  degree.  These  satisfy  the 
relations 


(A.  5) 

(A.6) 
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HK+ift)=  2^HKteWKHK_llZ>. 


TABLE  1 


Case 

Cl 

C2 

C3 

Clf 

Ql 

D2 

£>4 

I 

1 

1 

l 

l 

i 

1 

1 

1 

II 

1 

1 

l 

l 

1.5 

1.5 

1.5 

1.5 

III 

1 

1 

l 

l 

2 

2 

2 

2 

IV 

1 

1 

l 

l 

1 

1.5 

1.5 

2 

V 

1 

1 

l 

l 

1 

l 

2 

2 

VI 

1 

1 

l 

l 

1. 

1.5 

1-5 

1.5 

VII 

2 

1 

1 

l 

1 

1.5 

1.5 

1.5 

VIII 

5 

1 

1 

l 

1 

lf\ 

• 

rH 

ir\ 

• 

rH 

1.5 

IX 

1 

1 

l 

l 

2 

1 

1 

l 

X 

2 

1 

1  ' 

l 

2 

1 

1 

1 

XI 

5 

1 

1 

l 

2 

1 

1 

1 

Component  Error  Coefficients  (Cj)  and  Tol¬ 
erances  (Dj  in  the  Eleven  Cases  Studied 
J 
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IA3LE  2 


Case 

sr 

m2 

M4 

I 

2 

1.165 

3.709 

II 

2 

2.206 

13*53 

III 

2 

3.095 

27.22 

IV 

2 

2.168 

13.10 

V 

2 

2.130 

12.67 

VI 

2 

1.946 

10.47 

VII 

2.646 

2.819 

21.60 

VIII 

5-292 

8.933 

182.5 

IX 

2 

1.647 

7.490 

X 

2.646 

3.968 

40.90 

XI 

5-292 

20.22 

988.5 

The  Standard  Deviation,  Sj,  of  the  Output 
Distribution  Calculated  hy  the  Normal  Re¬ 
lation,  and  the  Output  Moments  Mg,  M^,  in 
the  Eleven  Cases  Studied 
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TABLE  3 


Case 

B0 

P* 

Monte 

Carlo 

P* 

Char. 

Function 

P* 

Mom. 

App. 

(I) 

1.6 

.8 

.822 

.858 

.858 

(III) 

2.8 

1.4 

.860 

.887 

.837 

(IV) 

2.4 

1.2 

.870 

.896 

.896 

(V) 

1.2 

.6 

.603 

•  579 

.580 

(V) 

2 

1.0 

.849 

.825 

.825 

(VI) 

1.2 

.6 

.612 

•  599 

.600 

(VIII) 

5 

.945 

.928 

.915 

.904 

(IX) 

1.2 

.6 

.677 

.639 

.640 

(IX) 

2.4 

1.2 

.912 

.940 

.940 

(X) 

3.175 

1.2 

.907 

.886 

.886 

Comparison  among  the  Values  of  P*  given 
hy  the  Three  Methods  of  CalcuJ.ation  for 

Various  Cases  and  Tolerances 
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0 


Figure  1 


2.0 


P*  (B0/Sj)  in  Case  I.  The  curve  is  ‘.he 
moment  approximation,  (^0),  and  the  circles 
are  Monte -Carlo  estimates 
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■5T— r~ 


Figure  6:  P*  (Bq/Sj)  in  Case  VI.  The  curve  is  the 

moment  approximation.  (40) ,  and  the  circles 
are  Monte -Carlo  estimates 

40 

* 

« 


Figure  7:  P*  (B0/Sj)  in  Case  VII.  'Tie  curve  is  the 

moment  approximation,  (40),  and  the  circles 
are  MonJ  ,  Carlo  estimates 


.5 


1.0 


1.5 

IVSj 


Figure  8:  P*  (B0/Sj)  in  Case  VIII.  The  curve  is  the 
moment  approximation,  (4o),  and  the  circles 
are  Monte -Carlo  estimates 
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Figure  10:  P*  (B0/Sj)  in  Case  X.  The  curve  is  the 

moment  approximation,  (40),  and  the  circles 
are  Monte -Carlo  estimates 
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